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ABSTRACT 

This paper describes the development of an LS-DYNA simulation of a vertical drop test 
of an ATR42-300 twin-turboprop high-wing commuter-class airplane. A 30-ft/s drop test 
of this aircraft was performed onto a concrete impact surface at the FAA Technical 
Center on July 30, 2003. The purpose of the test was to evaluate the structural response 
of a commuter-class aircraft when subjected to a severe, but survivable, impact. The 
aircraft was configured with crew and passenger seats, anthropomorphic test dummies, 
forward and aft luggage, instrumentation, and onboard data acquisition systems. The 
wings were filled with approximately 8,700 lb. of water to represent the fuel and the 
aircraft weighed a total of 33,200 lb. The model, which consisted of 57,643 nodes and 
62,979 elements, was developed from direct measurements of the airframe geometry, 
over a period of approximately 8 months. The seats, dummies, luggage, fuel, and other 
ballast were represented using concentrated masses. Comparisons were made of the 
structural deformation and failure behavior of the airframe, as well as selected 
acceleration time history responses. 

INTRODUCTION 

This paper describes the development of a full-scale finite element model of an ATR42 
commuter-class aircraft for crash simulation. The model was developed prior to the test 
and pre-test predictions were generated. Model validations are necessary to gain 
confidence in the application of explicit transient dynamic finite element codes for 
crashworthy design and certification. In fact, the “validation of numerical simulations” 
was identified as one of five key technology shortfalls during the Workshop on 
Computational Methods for Crashworthiness that was held at NASA Langley Research 
Center in 1992 [1], 

In 1998, the U.S. Army Research Laboratory Vehicle Technology Directorate (ARL- 
VTD) entered into an Inter-Agency Agreement (I A A) with the Federal Aviation 
Administration (FAA) William J. Hughes Technical Center for the purpose of validating 
crash simulations of airframe structures. As part of the IAA, finite element models were 
constructed of two 10-ft. long Boeing 737 (B737) fuselage sections, one outfitted with an 
auxiliary fuel tank mounted beneath the floor [2] and the other with two different 
overhead stowage bins and luggage [3]. Vertical drop tests of these two fuselage sections 
were performed at the FAA Technical Center in 1999 and 2000, respectively. These tests 
provided a valuable opportunity to evaluate the capabilities of computational tools for 



crash simulation through analytical/experimental correlation. Full-scale three- 
dimensional finite element models of the B737 fuselage sections were developed using 
MSC.Dytran [4], a commercial explicit transient dynamic code, and simulations of the 
vertical drop tests were executed. The analytical predictions were successfully validated 
through extensive test-analysis correlation, as documented in References 5 through 7. 

In 2003, the IAA was extended for an additional five years (through 2008) and the model 
validation work entered a new phase with the development of a full-scale finite element 
model of the ATR42-300 aircraft. For this simulation, the model was developed using 
the pre-processing software package, MSC.Patran [8], and the final model was executed 
using another commercial code, LS-DYNA [9]. The FAA performed a 30-ft/s vertical 
drop test of the aircraft to determine the impact response of the airframe; the floor, seat 
tracks, seats, and Anthropomorphic Test Dummies (ATD); and, the high-wing fuel 
system. Recently, the FAA has proposed dynamic performance criteria for seats for 
commuter-class aircraft that were based on empirical information obtained from prior 
airplane crash test data, which did not include airplanes representative in size of 
commuter-class aircraft. Consequently, this experiment was performed to provide impact 
data to evaluate the seat standards for commuter-class aircraft. The experimental 
program, the model development process, and the test-analysis correlation are presented 
in subsequent sections of the paper. 

EXPERIMENTAL PROGRAM 

On July 30, 2003, a 30-ft/s vertical drop test of an ATR42 aircraft was conducted using 
the Dynamic Drop Test Facility located at the FAA Technical Center. A pre-test 
photograph of the test article, raised to the drop height of 14-ft., is shown in Figure 1. 
This twin-turboprop, high-wing, commuter-class aircraft was developed and 
manufactured through a joint effort by Aerospatiale in France and Aeritalia in Italy. The 
aircraft has a wingspan of 80 ft., a seating capacity of 42-50 passengers, a maximum 
cruise speed of 304 knots/hour, and a maximum gross take-off weight of approximately 
41,000 lb. The drop test was performed onto a concrete surface. The purpose of the test 
was to evaluate the dynamic structural response of the aircraft when subjected to a 
severe, but survivable, impact. Particular attention was given to the seat and occupant 
responses to evaluate the FAA’s proposed dynamic seat requirements for commuter-class 
aircraft. 

A schematic drawing of the floor of the aircraft is shown in Figure 2. The total weight of 
the aircraft prior to the drop test was 33,200 lb. A large portion of the total weight was 
the 8,700-lb. of water added to the fuel tanks in the wings to represent fuel. In addition, 
16 double-occupant aircraft seats weighing 51 lb. each and 3 crew seats weighing 20 lb. 
each were attached to seat tracks on the floor, as shown in Figure 2. Seven ATDs and 16 
mannequins, each weighing approximately 170 lb., were seated in various locations, as 
indicated in Figure 2. Ballast weights equaling 135- and 152-lb. were added to some of 
the empty seats to represent the mass of occupants. The forward and aft storage 
compartments were filled with 1,450- and 739-lb. of luggage, respectively. In addition to 
the ballast, two overhead stowage bins were mounted between Frame Station 30 (FS 30) 



and FS 34 on the right and left sides of the fuselage. The bins weighed 55 lb. each (22-lb. 
empty weight plus 33-lb. ballast in each bin). Two simulated engine masses were added 
to the engine mounts on the wings, each weighing 1, 290 lb. 

A dominant feature of this aircraft is the high wing, which is attached directly to heavy 
fuselage frames at FS 25 and FS 27 through four dog bone specimens. The dog bone 
specimens attach to the fuselage frames at a location approximately 60 in. above the 
floor. Additional attachments are located at the very top of the fuselage frames. 

The measured longitudinal position of the center-of-gravity (CG) of the aircraft is 
approximately half way between FS 25 and 26, as shown in Figure 2. This measurement 
closely matches the CG location specified by the airframe manufacturer, as indicated in 
the Weight and Balance Manual [10]. 



Figure 1. Pre-test photograph of the ATR42 aircraft, raised to a drop height of 14 feet. 

The fuselage section was instrumented with accelerometers, strain gages, load cells, 
pressure transducers, and string potentiometers. Test data were collected at 10,000 
samples per second using an onboard digital data acquisition system. Of the channels 
available, eight acceleration responses were selected for correlation with the model. Of 
these eight responses, five were from accelerometers mounted on the floor of the fuselage 
cabin at locations shown in Figure 2. In addition, accelerometers located in the tail 
section at FS 47, the left sidewall at FS 18, and the center ceiling at FS 26 were also 
selected for comparison. These locations are not shown in Figure 2. 

A photograph showing an overall exterior view of the aircraft post- test is shown in Figure 
3. The dominant damage mode to the airframe was the failure of the heavy fuselage 
frames supporting the wing, causing the wing to subsequently displace through the 
fuselage cabin. The aluminum structure supporting the wing was crushed and fractured. 
A photograph showing a closer view of the wing/fuselage region post-test is shown in 
Figure 4(a). Further inspection of the airframe following the test indicated that several of 
the seats failed, as shown in Figure 4(b). 
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Figure 2. Floor schematic. 
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Figure 3. Photograph of the ATR42 aircraft post-test 





(a) Close-up view. (b) Interior view. 

Figure 4. Post-test photographs. 

MODEL DEVELOPMENT 

The finite element model of the ATR42 aircraft was developed from geometric data 
gathered from direct measurements of the aircraft, which were input into MSC.Patran 
database files. The geometric data were received from the FAA in six installments, each 
containing increasingly more detailed information on the shape and dimensions of the 
airframe. As each file was received, the new information was input into a master 
geometry file of the entire aircraft. The final geometry model of the aircraft, shown in 
Figure 5, consisted of 25,917 points; 17,270 curves; and 17,768 surfaces. The geometry 
model was discretized into a finite element mesh, element and material properties were 
assigned, contact and initial velocity conditions were defined, and the model was 
executed to generate analytical predictions of structural deformation and time history 
responses. 


Figure 5. Final geometry model of the ATR42 aircraft. 



The finite element model of the ATR42 aircraft, shown in Figure 6, contained 57,643 
nodes and 62,979 elements including 60,197 quadrilateral shell elements; 551 triangular 
shell elements; 526 beam elements; and, 1,705 point elements. 



(a) Front view. 



(b) Side view. 



(c) Top view. 



(d) Three-quarter view. 


Figure 6. Finite element model of the ATR42 aircraft. 


An automatic contact (CONTACT AUTOMATIC SINGLE SURFACE) was specified 
for the model, which is a generic contact definition that prescribes that no node can 
penetrate through any surface in the model. An impact surface was created to represent 
the concrete pad beneath the drop tower. This surface was modeled as a 5-in. thick 
aluminum plate, encompassing the total length and width of the aircraft, as shown in 



Figure 6(d). Three main material properties are defined in the model for aluminum 
A12024-T3, aluminum A17076-T5, and titanium Ti-6A1-4V. These properties were 
defined using the MAT PLASTIC KINEMATIC card in LS-DYNA for a linear elastic- 
plastic material with input values for density, Poisson’s ratio, Young’s modulus, yield 
stress, hardening modulus, and an ultimate failure strain. Most of the sheet metal parts, 
such as the outer skin were assigned material properties of A12024-T3. The forged metal 
parts, such as the fuselage frames, floor beams, and seat tracks were assigned material 
properties of A17076-T5. The “dog bone” specimens used to attach the wing to the 
fuselage frames at FS 25 and FS 27 were assigned material properties of titanium. The 
material property designations for each component were obtained from the aircraft 
manufacturer’s Weight and Balance Manual [10] and the material property values were 
found in MIL-HDBK-5H [11]. 

The pictures of the model, shown in Figure 6, were obtained from FS-POST [12], the 
post-processing software for FS-DYNA. A figure depicting the location of the point 
elements is shown in Figure 7. Point elements were used to assign concentrated masses 
representing the seats, occupants, luggage, fuel, and other ballast to nodes in the model. 
The mass distribution of the fuel loading was determined by creating solid elements 
defined by the upper and lower wing skins, with one element thickness between each 
skin. The nodes on the bottom skin were fixed. The solid elements were assigned 
properties with the density of water, and a gravitational loading was applied. The 
resulting constraint forces on the bottom nodes were converted into a mass loading at 
each node. Once this process was complete, the solid elements and nodal constraints 
were then removed. 


All of the nodes in the aircraft model were given an initial velocity of -360 in/s (-30 ft/s). 
For the impact surface, all of the edge nodes were constrained from translational and 
rotational motion. The model was executed in FS-DYNA (version 970) for 0.3 seconds 
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Figure 7. The location of point elements in the model. 


of simulation time, which required 130 hours of CPU on a single processor Hewlett 
Packard workstation x4000. 

As a quality check on the model, the total weight and longitudinal CG location of the 
model were compared with those of the test article. The weight of the aircraft was 33,200 
lb. and the total weight of the aircraft model was 32,900 lb., just 300 lb. lighter than the 
test article. The measured longitudinal CG position of the test article was 469.2 in. from 
the reference location. For the model, the longitudinal CG position was 474. 1 in. from 
the reference location, within 5 inches of the experimental value. 

TEST-ANALYSIS CORRELATION 

The test-analysis correlation consists of comparisons of experimental and predicted 
structural deformation and selected acceleration and velocity time history responses. 

Comparison of Structural Deformation 

A comparison of structural deformation is provided in Figure 8, which shows a 
comparison of test article and model deformation from 0.06 to 0.3 seconds in .06-second 
intervals. 



Figure 8. Comparison of structural deformation. 


The pictures of model deformation were obtained from the post-processing file, and the 
experimental pictures were captured from the high-speed video. In general, the model 
accurately predicts the structural deformation and failure behavior of the test article, 
including collapse and failure of the fuselage structure beneath the wing. In both the test 
and analysis, the failure is initiated by fracture of the fuselage frames at FS 25 and FS 27, 
and not failure of the dog bone specimens. The frame failures allow the wing to translate 
downward through the fuselage cabin. However, some differences are noted. In the 
model, the wing exhibits considerable flexure and bending at the tips and this behavior is 
not seen in the test. 

Comparison of Selected Acceleration Time History Responses 

The raw experimental acceleration data were plotted versus time, integrated to obtain the 
velocity time history response, and filtered using a low-pass digital filter based on the 
SAE J211/1 specifications [13]. A variety of filter frequencies were applied in an 
attempt to obtain a discernable single acceleration pulse, while at the same time not 
distorting the integrated velocity response. The distortion is determined by comparing 
the velocity responses obtained by integrating the raw acceleration data and by 
integrating the filtered acceleration data. A filtering frequency of 33.2 Hz was selected 
because it met these criteria. The experimental and analytical acceleration pulses shown 
in this section of the paper were filtered at this frequency. 

The experimental and analytical acceleration responses for the left and right outer seat 
track locations at FS 20 are shown in Figure 9. For the test, the main acceleration pulse is 
less than 0.1 -second in duration, and exhibits two main peaks. For the left outer seat 
track location, the first acceleration peak is smaller in magnitude than the second, which 
has a magnitude of 28 g. For the right outer seat track location, the two acceleration 
peaks are of the same magnitude, about 24 g. In both cases, the predicted acceleration 
responses also exhibit two peaks, with the second being higher in magnitude than the 
first. The pulse durations of the analytical responses are almost identical to that of the 
experimental responses. For comparison, the raw experimental and analytical 
acceleration data were integrated to obtain the velocity time history responses at these 
two locations. The velocity responses are shown in Figure 10. The level of agreement 
between the test and analysis is good, though the test data shows a higher rebound 
velocity than predicted by the analysis for the right outer seat track location. 

Next, the experimental and analytical acceleration responses of the right side of the 
cockpit floor are plotted in Figure 11. Again, the experimental acceleration response 
exhibits two peaks; however, in this case the first peak of 34 g is higher in magnitude 
than the second of 28 g. The predicted acceleration response shows a single pulse, of 
longer duration than the experiment, with a peak acceleration of 32 g. 

The experimental and analytical acceleration responses of the left and right inner seat 
track locations at FS 29 and FS 35, respectively, are plotted in Figure 12. The locations 
of these two accelerometers are shown in the schematic drawing of Figure 2. The 
accelerometer at FS 29 is located slightly to the rear of the fuselage frames supporting the 
wings, while the accelerometer at FS 35 is located at the very rear of the aircraft. The 



filtered experimental accelerations exhibit high-amplitude, low frequency responses, 
making it difficult to discern a single acceleration pulse. In general, the predicted 
acceleration responses show good agreement with the experimental data at these two 
locations. 
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(a) Left outer seat track at FS 20. (b) Right outer seat track at FS 20. 

Figure 9. Experimental and analytical acceleration responses for the left and right outer 

seat track locations at FS 20. 
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(a) Left outer seat track at FS 20. (b) Right outer seat track at FS 20. 

Figure 10. Experimental and analytical velocity responses for the left and right outer 

seat track locations at FS 20. 
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Figure 11. Experimental and analytical acceleration responses of the right side of the 

cockpit floor. 


The experimental and predicted acceleration responses of the center of the tail section at 
FS 47 are plotted in Figure 13. Unlike the floor acceleration responses, which had pulse 
durations of less than 0.1 second, this response is 0.2 seconds long. Both the 
experimental and analytical acceleration responses exhibit a single pulse, of 
approximately the same duration, with magnitudes of 12 and 13 g’s, respectively. 
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(b) Right inner seat track at FS 35. 


Figure 12. Experimental and analytical acceleration responses of the left and right inner 
seat track locations at FS 29 and FS 35, respectively. 
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Figure 13. Experimental and analytical acceleration responses of the center of the 

fuselage at FS 47. 

The experimental and analytical acceleration responses of the left sidewall at FS 18 are 
plotted in Figure 14. This accelerometer was located on the sidewall approximately 12 
in. above the floor and was oriented in the vertical direction. The analytical acceleration 
response closely matches the magnitude (peak acceleration of 25 g’s for the analysis 
compared with 22.5 g’s for the experiment) and duration of the experimental pulse. 
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Figure 14. Test and analysis acceleration responses of the left sidewall at FS 18. 

The final time history comparison is shown in Figure 15, in which the experimental and 
analytical acceleration responses of the center ceiling at FS 26 are plotted. This location 
is essentially at the center of the bottom skin of the wing. The filtered experimental and 




analytical data continue to exhibit high-amplitude, low frequency responses at this 
location and it is not possible to discern a single acceleration pulse. The predicted 
response overshoots the magnitude of the initial peak acceleration; however, it accurately 
captures the dip and subsequent rise in the experimental response that occurs at 
approximately 0.1 seconds. 

Acceleration, g 
20 

15 
10 
5 
0 
-5 
-10 

- 0.05 0 0.05 0.1 0.15 0.2 0.25 

Time, s 

Figure 15. Test and analysis acceleration responses of the center ceiling at FS 26. 

DISCUSSION OF RESULTS 

In general, a high level of agreement was obtained between the experimental and 
analytical data, especially when considering the complexity of the test article and the 
unsophisticated method of model development. The simulation accurately predicted the 
major structural failure, the collapse and failure of the fuselage frames supporting the 
wing structure. However, some differences in the progression of structural damage and 
deformation were noted. The model exhibited more bending and flexure of the wing than 
was observed experimentally, indicating that the material properties used in the model 
were inaccurate. 

The model predicted the experimental acceleration time histories quite well. The good 
agreement achieved for the floor-level acceleration responses is important in that these 
pulses are subsequently transmitted to the seats and occupants and this information is 
required to establish dynamic seat criteria for computer-class aircraft. It was noted that 
even after filtering, some of the acceleration responses exhibited high- amplitude, low 
frequency responses, making it difficult to discern a single acceleration pulse. For other 
channels, the low-pass filter with a cut-off frequency of 33.2 Hz was able to generate a 
single acceleration pulse. Thus, for this test, it may be necessary to use different cut-off 
frequencies, on a channel-by-channel basis. Of course, the analytical data would be 
filtered at the same frequency as the experimental data at each location. 





CONCLUDING REMARKS 


A full-scale three-dimensional finite element model of a twin-turboprop high-wing 
commuter-class aircraft, the ATR42-300, was developed and executed as a crash 
simulation. The analytical predictions were correlated with test data obtained from a 30- 
ft/s vertical drop test of the aircraft that was conducted using the Dynamics Drop Test 
Facility at the FAA William J. Hughes Technical Center in Atlantic City, NJ. For the 
test, the aircraft was configured with seats, anthropomorphic test devices, luggage in the 
forward and aft compartments, and 8,700-lb of water in the wings to represent the fuel 
loading. The finite element model of the aircraft was developed from direct geometric 
measurements and contained 57,643 nodes and 62,979 elements including 60,197 
quadrilateral shell elements; 551 triangular shell elements; 526 beam elements; and, 
1,705 point elements. The model was executed in LS-DYNA, a commercial code for 
performing explicit transient dynamic simulations. 

The analytical predictions correctly simulated the major damage mode seen during the 
test which was collapse and failure of the fuselage structure beneath the wing. These 
structural failures allowed the wing to displace vertically through the fuselage cabin. 
Finally, the predicted acceleration responses showed a high level of correlation with the 
test data, which is remarkable considering the complexity of the test article and the 
unsophisticated method of model development. It was particularly important to obtain 
accurate prediction of the floor-level acceleration responses, since these pulses are 
transmitted to the seat and occupants. This data will be useful in evaluating the FAA’s 
proposed dynamic seat standards for commuter-class aircraft. 
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